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Abstract 



Using general features of recent quantizations of the Hamiltonian constraint in 
loop quantum gravity and loop quantum cosmology, a dynamical interpretation of 
J> ' the constraint equation as evolution equation is presented. This involves a trans- 

■ formation from the connection to a dreibein representation and the selection of an 

'. internal time variable. Due to the discrete nature of geometrical quantities in loop 

QO quantum gravity also time turns out to be discrete leading to a difference rather 

than differential evolution equation. Furthermore, evolving observables are discussed 
^3 ■ in this framework which enables an investigation of physical spectra of geometrical 

Q I quantities. In particular, the physical volume spectrum is proven to equal the discrete 

O^' kinematical volume spectrum in loop quantum cosmology. 

U 

bp. 

1 Introduction 

X 

i-j I One of the major open issues of loop quantum gravity |l|] is to understand its dynamics 
which is governed by the Wheeler-DeWitt equation. However, this equation in the full 
theory (e.g. 0) contains all the complicated details of the evolution and is plagued not 
only by technical problems but also by conceptual issues, most importantly the problem of 
time. At this point, reductions to simpler but still representative models may be helpful, 
a strategy which in previous approaches to (quantum) gravity has been widely used. Dy- 
namical issues are most conveniently studied in cosmological models which are obtained by 
a reduction of gravity to homogeneous geometries. For loop quantum gravity, this reduc- 
tion has been performed in Ref. by specializing the general symmetry reduction scheme 
for diffeomorphism invariant quantum theories of connections [Q. In this framework we 
already derived quantizations of the Hamiltonian constraint for various models which 
will be used now to investigate the dynamics governed by those operators. 

*e-niail address: bojowald@physik.rwth-aachen.de 
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The study of (quantum) cosmological models has been initiated in order to understand 
very early stages of the universe and to gain insights which can also be used in the full 
theory. In lack of a complete quantum theory of gravity, however, the only route to quan- 
tum cosmology has been to perform a quantization after the classical symmetry reduction 
(minisuperspace quantization @, Although it is not at all clear, and can from this 
perspective not be decided, whether the immense restriction to finitely many degrees of 
freedom is mild enough to maintain typical properties of the full quantum gravity, minisu- 
perspace models provide interesting test-beds for problems of quantum gravity. Due to the 
finite number of their degrees of freedom they are quantum mechanical models, but their 
dynamics is, as inherited from the theory of General Relativity, still intrinsic. This means 
that it is not governed by a Schrodinger equation for the evolution of a wave function in 
an external time parameter, but by the Wheeler-DeWitt equation which is a constraint 
equation solely expressed in terms of the metrical variables and their conjugate momenta 
corresponding to the fact that the theory is invariant under arbitrary reparameterizations 
of time. At first sight, there seems to be no dynamics at all in this constraint formulation, 
which lead to the name "frozen time formalism" . 

In order to interpret the Wheeler-DeWitt equation as a time evolution equation one 
has to introduce an internal time which is constructed from metrical or matter degrees of 
freedom (see, e.g., Ref. and references therein). In a homogeneous context one usually 
chooses the scale factor of the universe (related to the determinant of the metric), in which 
the Wheeler-DeWitt equation is a hyperbolic differential equation of second order (this 
holds true also after introducing inhomogeneous perturbations 0). 

Let us illustrate these considerations with simple examples. In an isotropic flat model, in 
which space looks the same not only in all points but also in all directions, the only metrical 
degree of freedom is described by the scale factor a appearing in the spatially isotropic 
metric ds^ = —dt^ + a(t)^(da;^ + dy^ + dz"^). A model with only one degree of freedom can, 
in a formalism without extrinsic time, not exhibit any dynamics: there simply is nothing 
to build a clock. In the space-time picture this corresponds to the fact that the classical 
solutions are maximally symmetric (not only in space but in space-time) : Minkowski space 
for a vanishing cosmological constant or DeSitter/Anti-DeSitter space. The most simple 
matter field which can be coupled is an isotropic scalar field 0. We now have two degrees 
of freedom, a and (p, and dynamics which is usually written down in terms of Lagrangian 
or Hamiltonian equations of motion. Although these equations are differential equations 
for a{t) and (p{t) in terms of a time parameter t, this parameter can be reparameterized 
(gauged) arbitrarily. In this simple model it is immediate to remove t in order to arrive at 
an intrinsic time formalism: from the differential equations for a{t) and (j){t) one can obtain 
a differential equation for ^(a), where a is regarded as intrinsic time. Its solution describes 
the evolution of the scalar field in an expanding or contracting branch of the universe 
and contains all invariant information about the model. In a more complicated situation 
one can, similarly, make sense only of relational motions of degrees of freedom with respect 
to each other, and not with respect to an external time. Analogously, in a minisuperspace 
quantization we can describe the states by wave functions ip{a, (j)) (often dubbed the "wave 
function of the Universe") subject to the Wheeler-DeWitt equation = 0. Here, no 
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external time parameter appears from the outset, and the information contained in -0(0, (p) 
is relationah interpreted as a probabihty density, it describes the possible values of 4> in 
relation to the values of a. 

Other widely used cosmological models are the Bianchi models, which describe homo- 
geneous, but not necessarily isotropic, geometries. Therefore, the metric is parameterized 
by more than one degree of freedom, and one can study dynamics of the vacuum solutions 
without coupling matter fields. This has been done in Ref. [|10| with the following basic 



results: The metric can be consistently diagonalized (in some models this is just a gauge 
fixing, whereas in others it is a further truncation in addition to the symmetry requirement 
111 ) such that there are only three coordinates usually denoted as /3° (the scale factor) 



and P^, (the anisotropics) and their canonically conjugate momenta ttq, tti, 7^2- In a 
suitable gauge which also specifies the lapse function the Hamiltonian is given by 

n = Iv'^nmj + V{f3') ^ 

where t] is a constant metric on minisuperspace with signature (— , +, +) and V a potential 
which characterizes the specific Bianchi model. For Bianchi I the potential V vanishes, 
which can also be achieved for other models after suitable coordinate transformations on 
minisuperspace, and the Hamiltonian consists of just the "kinetic" term containing the 
momenta. 

In the quantum theory of this model we can choose the /3-representation in which wave 
functions ip depend on the parameters , i.e. they are functions on minisuperspace. As 
usual, the momenta then are represented as derivative operators and the Wheeler-DeWitt 
equation takes the form 

of a Klein-Gordon equation where P^, again the scale factor, plays the role of time. An 
interpretation as evolution equation in the intrinsic time is then immediate. However, 
we can just as well choose a 7r-representation by using wave functions which depend on the 
momenta. In this case, the Wheeler-DeWitt operator would not be a derivative but a mul- 
tiplication operator, and the Wheeler-DeWitt equation would not have an interpretation 
as evolution equation but instead constrain the support of wave functions. Of course, both 
pictures are equivalent, as in usual quantum mechanics, but we see that the emergence 
of an evolution equation depends on the representation once an internal time is selected. 
This issue is characteristic for a generally covariant model where the internal time is to be 
found under the internal degrees of freedom, which do not have a unique representation. 

There is an important lesson we have to learn from these considerations: In cosmological 
scenarios it is most convenient to choose the scale factor as time variable (although this 
was the case in both examples, the selection of an internal time is by no means unique; 
furthermore, in more complicated models, let alone the full theory, there are no explicitly 
known time functionals) and to study the evolution of other degrees of freedom (matter or 
metrical) with respect to this parameter. In a minisuperspace quantization one then has 
to use a metric representation in order to extract an evolution equation. Using metrical 
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variables {qab^p""^) this is the usual representation anyway, but the largest successes with 
respect to kinematical aspects of quantization have been achieved in loop quantum gravity 
using connection variables where one bases the quantum theory on the connection (or loop) 
representation. Also the quantum symmetry reduction to homogeneous models leads at 
first to a connection representation such that, afterwards, we have to transform to a dreibein 
representation and find an interpretation of the Wheeler-DeWitt equation as evolution 
equation^. Here, one can expect significant departures from the usual minisuperspace 
quantizations described above, because the volume, which has been used as internal time. 



is now quantized at least at the kinematical level [|T2], Q (note that the parameters a and 



/3° above can take all positive values and thus show no volume quantization). We will see 
that this implies a discrete time and a difference (not differential) equation as evolution 
equation. 

Our strategy is the following: Using Ref. |^ we perform the symmetry reduction 
for the above models at the kinematical level of loop quantum gravity by selecting homo- 
geneous states. Alternatively, the procedure can be interpreted as a loop quantization of 
minisuperspaces. As shown in Ref. |jl5[, this leads to discrete geometric spectra after the 



symmetry reduction, whereas a usual minisuperspace quantization after a classical sym- 
metry reduction leads to a continuous volume spectrum. The dynamics of our models is 
governed by the Hamiltonian constraint operators which have been derived in Ref. and 
will be recalled briefly in the following section. In the main part of this paper we will 
perform the transformation into a dreibein representation and discuss our interpretation 
of the constraint equation as an evolution equation. Finally, having an evolution equa- 
tion, we can study evolving observables and will, in particular, show that the physical 
volume spectrum (taking into account the Hamiltonian constraint) equals the kinematical 
one (ignoring the Hamiltonian constraint). 



2 Wheeler— DeWitt Operators 

We will first recall the loop quantum theory of homogeneous minisuperspaces derived in 
Ref. and the Hamiltonian constraint operators of those models [^] which will be used 
in the remaining part of this article. 

Because of their homogeneity the values of all fields in a single point suffice to completely 
characterize the spatially homogeneous canonical fields, namely the scalars 0}ri G CSU{2), 
1 < / < 3 which determine a homogeneous connection Al^Ti = (p^TiUO^ in terms of invariant 
(with respect to the symmetry group) one-forms on E and their canonically conjugate 
momenta which are derived from the dreibein components and, therefore, encode the 
degrees of freedom of the metric on E. Consequently, a reduced formulation can be formu- 
lated in a single point rather than by fields on a space manifold S which explicitly shows 
the reduction to finitely many degrees of freedom. This also demonstrates that usual loop 
variables cannot be used to build an auxiliary Hilbert space because they would necessarily 

^The author is grateful to A. Ashtekar for this suggestion. 
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break the symmetry. Instead, one uses point holonomies which are SU (2)-elements asso- 
ciated with a single point and serve to describe quantum scalar fields. For Bianchi models 
(because we will use a Hamiltonian formulation we have to restrict our considerations to 
class A models whose structure constants cfj fulfill cjj = 0) there are three independent 
point holonomies hj G SU{2), 1 < / < 3 (due to anisotropy) and quantum states are 
gauge invariant functions (where the S'[/(2)-gauge group acts by joint conjugation of all 
three point holonomies) of these three group elements. The behavior under conjugation 
shows that a point holonomy can be represented as an ordinary holonomy associated to a 
closed loop embedded in an auxiliary manifold, and quantum states can be expanded in 
spin network states on graphs with three closed edges e/ meeting in a 6-vertex. We are 
then able to compare results derived for minisuperspaces with those of the full theory by 
restricting the latter to spin networks with a single 6-vertex. For instance, the volume 



operator for Bianchi models is identical to the action on a 6-vertex of the volume 
operator in the full theory derived in Ref . [T^ (but differs from that of Ref . |]12| , see also 
Ref. m). 



For locally rotationally symmetric (LRS) or isotropic models, the three point holonomies 
are no longer independent, but related by the equations 

Mp(/)(e,))=Ad,(^)Me/) (1) 

where / is an element of the isotropy group acting on the edges ei via the representation 
p, and A is a homomorphism from the isotropy group to the gauge group SU (2) (see Ref. 

for details). In particular, for isotropic models all three point holonomies are related 
and quantum states can be expressed as functions on a single copy of SU{2). However, 
they are not ordinary spin network states on a single edge, but generalized spin network 
states which can have an insertion due to the gauge transformation on the right hand side 
of Eq. (0). Details have been worked out in Ref. [jl5l where also the volume operator has 
been derived and diagonalized explicitly. 

Finally, we will need the Hamiltonian constraint operators which can be derived 
similarly to the full theory [^. Adaptations in the regularization are necessary only be- 
cause there is no continuum limit and because the reduced operators have to respect the 
symmetry. But the splitting of the Lorentzian constraint in a Euclidean part and a poten- 
tial term, and the usage of the extrinsic curvature can be adopted without changes. For 
Bianchi models the Euclidean constraint operator is 

n^^'^N] = -Ai{i'll)-^NY,^''^^^ {hihjh-,%'h-^,]j^[hKM) (2) 

UK 

with the Planck length l-p and l' = lVq'^ being related to the Immirzi parameter l and the 
volume Vq of space in a fiducial homogeneous metric. Furthermore, hj are point holonomies 
acting as multiplication operators using the definition 



K=l 
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and V is the volume operator for Bianchi models. 
Using the quantized extrinsic curvature 



k = itr^ 



the complete Lorentzian constraint can be written as 

7t:[A^] = 8^(1 + ^2) (4)-3yoA^e^"^'^trf[/i,,K][/ij,i^][/i;^,y]) -U'^^^N] 



(3) 



Constraint operators for LRS and isotropic models can be derived by inserting the 
conditions ([l|) into the operators for Bianchi models and evaluating the action on the 
reduced Hilbert spaces. 



3 Dynamics 

Following the usual procedure for interpreting the Wheeler-DeWitt equation of cosmologi- 
cal models as an evolution equation, we will now first transform to a dreibein representation 
and then introduce an internal time. This will allow us to demonstrate that the Hamilto- 
nian constraint equation can be written as an evolution equation with discrete time. 

3.1 Dreibein Representation 

In most discussions of loop quantum gravity one works always in the connection represen- 
tation where quantum states are represented as functions on the space of connections which 
are usually expanded as linear combinations of spin network states. Also the Hamiltonian 
constraint has been quantized to an operator acting on this space of functions where the 
holonomies, which constitute the main part of the constraint operator, act as multiplication 
operators 0, ^ . Thus imposing the Wheeler-DeWitt equation will result in a restriction of 
the support of physical states being in the kernel of the quantum constraint. This is similar 
to the discussion recalled in the Introduction of the standard minisuperspace quantization 
of the Bianchi I model, where the Wheeler-DeWitt equation in the vr-representation re- 
stricts the support of wave functions. Note that the variables tt/ there are conjugate to the 
metrical variables which also are the connection variables in the present framework. As 
suggested in this context, formulated in a metric (or dreibein) representation the Wheeler- 
DeWitt equation usually allows an interpretation as evolution equation which motivates 
the following discussion. 

In a dreibein representation quantum states are represented as functions on the space of 
dreibein components. Whereas for the homogeneous models which are of interest here this 
space is finite-dimensional, in the full theory it is infinite-dimensional and a mathematical 
formulation on the relevant function spaces has to be done with great care. Up to now 
we always took the attitude that all our constructions in the reduced models should be as 
close to the full theory as possible, and using a technique which cannot be generalized to 
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the infinite-dimensional space of tlie full theory would obviously spoil this aim. A rigorous 
formulation of the connection representation has been achieved by applying the theory of 
representations of C*-algebras to a particular C*-algebra constructed from holonomies. A 
similar procedure is not applicable for a dreibein representation and, whereas the space of 
connections could be compactified to the quantum configuration space A of the connection 
representation, a suitable compactification of the space of dreibein components is not 
obvious and, even worse, not reasonable from the physical point of view. Furthermore, 
the compactness of the space A means that in a quantum theory the conjugate momenta, 
i.e. the dreibein components, can have only discrete values, which we have already seen 
in the discrete spectra of geometrical observables (the dreibein components are quantized 
to angular momentum operators). Thus the quantum configuration space in a dreibein 
representation will be a discrete space. 

We will follow here a strategy which makes use of the already known and mathematically 
well-founded connection representation by constructing a transform from this representa- 
tion to another representation which will be equivalent to a dreibein representation. As is 
well-known from quantum mechanics, such a transform can be constructed by expanding 
a state in a given representation in terms of eigenstates of a complete set of commuting 
operators. In our case, these operators should be quantizations of the dreibein components 
or, in order to maintain gauge invariance, of the products pfjo"^* which correspond to the 
metric components g^'^ . More convenient a procedure is to expand a state in the connection 
representation in spin network states (which are usually chosen as a basis, anyway) and 
to use the spins and vertex contractors (which can also represented by spins) as discrete 
coordinates of the "quantum dreibein space". Obviously, all the spins (in the full theory 
we would have to include other discrete labels, e.g. knot invariants, parameterizing the 
diffeomorphism equivalence classes of graphs) form a complete set characterizing a state 
completely, and this description is equivalent to a dreibein or metric representation because 
all eigenvalues of metric components can be expressed in terms of the spins, and vice versa. 

Introducing a model dependent index set X which contains all allowed multi-labels L, 
we can write the decomposition of an arbitrary state / in the connection representation as 

Lei 

where {Ti : L G X} is an orthonormalized set of spin network states associated to the 
multi-labels L. For the Bianchi models we have 

2^Bianchi = {L= (il,i2,i3, ^1, ^2, h) ■ 3U kj E ^Nq} 

if we parameterize the 6- vertex contractor by the three spins kj {ji are spins associated to 
the external edges of the point holonomies). For LRS models we have two external edges, 
a contractor parameterized by a single spin k, and an additional label i describing the 
insertion; for isotropic models there is only one spin j and an insertion-label i taking only 
two possible values (see Ref. p!5[] ): 

^LRS = {{3i,j2,k,i)} , Xi3o = {(j,0 : J e iNoU{-i},«G {0,1}} . 
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From now on we will call a state / represented by the coefficients Jl in the above 
expansion a state in the dreibein representation. In this representation states are maps 
/: X ^ C, L 1-^ /i from the index set of the respective model to the complex numbers. The 
inner product can be derived from that in the connection representation (the Ashtekar- 
Lewandowski inner product): because the spin network states in the expansion of f{A) 
were assumed to be orthonormalized in the Ashtekar-Lewandowski measure, we have 

if, 9)= h9L'iTL,TL>) = J2h9L. (5) 
L,L'ei Lex 

Thus, the kinematical Hilbert space is represented in the dreibein representation as Tikm = 
£^(X), the completion of the space of square integrable sequences on the index set X. 



3.2 Internal Time 

A central ingredient for a dynamical interpretation of a generally covariant theory is to 
choose one combination of the degrees of freedom as an internal time. In standard quanti- 
zations of homogeneous minisuperspaces (and their classical counterparts) this is usually 
done by using the scale factor (related to the spatially constant determinant of the met- 
ric on E) of the universe 0. Besides providing the intuition that time is related to the 
expansion (or contraction) of the universe, this has the virtue of resulting in a hyperbolic 
differential equation (Wheeler-DeWitt equation) governing the evolution by means of a 
well-posed initial value problem. 

In principle, we could copy this procedure and try to extract some interpretation of 
time from the volume spectrum in our quantizations of homogeneous models. However, 
the volume spectrum is, in general, quite complicated and not known explicitly for the 
Bianchi and LRS models. Although we know the volume spectrum for isotropic models 



15| , these models have only a single gravitational degree of freedom and we have to couple 
matter in order to obtain a reasonable dynamical system. 

We therefore look for an acceptable substitute of the volume as internal time which 
we will motivate by using the classical Bianchi I model on M^. Its solutions, the Kasner 



solutions ^7j, can be written as 

ds^ = -dt^ + f^'dxl + t^^^dxl + f'^^dxl 

where the parameters a, have to obey the relations X]/'^/ — X]/'^/ — i-^- there is 
only one independent parameter. One can see that there are always two positive and one 
negative parameter, so that two directions of space are expanding and the third one is 
contracting in such a way that the volume increases monotonically (due to Yli aj = 1). In 
other models this behavior is also the generic one in certain time intervals, the so-called 
Kasner epochs, where, however, a transition between different epochs is possible. The 
Kasner behavior demonstrates that we can expect each diagonal metric component g^^ , 
describing the expansion or contraction of the I-th direction, to be as good an internal 
time as the scale factor. In a standard minisuperspace quantization we would now have to 
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show that we get a hyperbohc evolution equation with such an internal time, but here we 
are at first mainly interested in a simple spectrum of the time parameter. 

Let us therefore pick the first metric component p\p^^ as a tentative time for Bianchi 
models. It is readily quantized using the usual procedure for point holonomies |^ : 

In order to determine its spectrum, we have to choose a suitable parameterization of the 
contractor in the 6-vertex, which can be done by decomposing it into a combination of 
four 3-vertices the contractors of which are unique up to a constant factor. We do this 
by contracting first each of the closed external edges carrying spins ji to an internal edge 
carrying spin ki\ the internal edges are then contracted in a central 3- vertex. We orient the 
edges in such a way that the internal ones are incoming in the central vertex, whereas the 
external edges have an incoming and an outgoing part in the three non-central vertices. 
All four 3-vertices are gauge invariant so that we have the relations 

Jf^\hj) - Jf'\hj) = Lf\l) 

defining Lf"\l) as a left invariant angular momentum operator associated with the J-th 
internal edge. 

The computation of the spectrum of the diagonal metric components is now similar to 
that of the area spectrum (and to computations of spectra of coupled angular momenta): 

PlP^^ = (2jW(/ii)2 + 2jW(/ii)' - L^'^\l?) 

immediately leads to the eigenvalues 

We now have a candidate for an internal time with an explicitly known and simple spec- 
trum, but we will even simplify this by using the external spin ji as time label. Although 
we did not justify it as labeling eigenvalues of a quantization of a classically admissible time 
function, it is favorable because its lowest value, ji = 0, implies vanishing volume V = Q 
(for a vanishing spin on one external edge the volume eigenvalues are those of a planar 4- 
vertex which always vanish: there are only two independent angular momentum operators 
from which no non- vanishing antisymmetric product in three indices can be built). Ulti- 
mately, the justification of an internal time has to come from a reasonable interpretation 
of the dynamics as evolution in that degree of freedom, which will be studied below. 

In other models, LRS or isotropic, we can similarly pick one of the labels of the quantum 
states as a candidate time label. We will in general decompose the label as (n, L) where n 
is the time label and L denotes labels for all other remaining metric or matter degrees of 
freedom. E.g. in Bianchi models, we have n = ji and L = {j2; is? ^i? ^2? ^s}? and in isotropic 
models n = j and L solely contains the insertion and possibly matter labels. Quantum 
states then are given in the dreibein representation by the coefficients Cn,L, which is a 
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discrete substitute of the wave function ip{a, 0) of standard minisuperspace quantizations. 
The decomposition of the labels corresponds to a decomposition of the kinematical Hilbert 
space 

HMn = f{I)=^Vn (6) 

n 

into "equal time" subspaces P„ in which the time label n is fixed, but the remaining labels 
in L are arbitrary. 

3.3 Discrete Time Evolution 

After transforming to the dreibein representation and picking a candidate for an internal 
time, we are now ready to study the dynamics governed by the Wheeler-DeWitt operator 
(0). The essential part of this operator is the multiplication with a couple of holonomies, 
where also the main contribution to the model dependence of the operator enters. 

In the connection representation the action of matrix elements of a holonomy as 
multiplication operator is given by 

where vr-' denotes the matrix representation of SU{2) associated with spin j. This action 
can be transformed to the dreibein representation and schematically written as (acting on 
a state c G f{I)) 

suppressing all but the one label associated with the edge underlying the holonomy h. 
Here, the coefficients hj are not just real numbers but operators acting on the subspaces 
of the kinematical Hilbert space with fixed j, i.e. in general they change the remaining 
labels which have been suppressed in the above equation. E.g. in the Bianchi models, the 
external spins ji, and therefore our internal time n = ji, are changed only by multiplication 
with the holonomy hj associated to the same edge as the spin. However, the internal spins 
kj are affected by all holonomy multiplications because they parameterize the contractor 
which always is subject to change. Similarly in LRS and isotropic models, external spins 
are changed only by multiplication with the appropriate holonomy, whereas the insertion 
and the contractor are changed by all holonomies. 

From now on we will mainly be interested in the "time spin" n and the holonomy 
(denoted as hn in what follows) affecting it. Note that this is possible only because we 
selected a spin of an external edge as our internal time; otherwise the following consid- 
erations would be more complicated. Each /i„ appearing in the constraint operator leads 
to a combination of coefficients with labels n + ^ and n — | of a state in the dreibein 
representation. Because there is always more than one holonomy associated with a fixed 
edge in the constraint operator, we have in general an action of the form 

(He) 
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again suppressing the labels in L and introducing the operators Hi which fix each of the 
subspaces Vn. The number cu is defined as twice the maximal number of holonomies hn 
appearing in each summand of the constraint operator; it also determines the number of 
operators Hi which is given by a; + 1. 

Because H is a quantized constraint operator, all physical states s have to obey the 
Wheeler-DeWitt equation 

2 

J2iH^sUi = ^ neiNo. (9) 

* 2 

In this form, recalling that n is our internal time label, we can immediately read off a dis- 
crete evolution equation: the Wheeler-DeWitt equation is an implicit difference equation 
of generic order ou which, however, is reduced up to ^ for small n because in the holonomy 
multiplication there is only one term when acting on the n = 0-state. Provided the highest 
order operator //| is invertiblc on each subspace Vn, wc can turn this into a difference 
equation where all coefficients c„,l are determined by the coefficients on the preceding to 
equal time subspaces, which directly leads to a formulation as an initial value problem 
which is uniquely specified by fixing initial conditions on the first ^ equal time subspaces 
Dq to Vu^_i (although the generic order is uj, for the initial value problem the order ^ for 
small n matters; note that n is a half-integer). If H^ is not invertible on some subspace 
Vn, then not all coefficients at time n are specified uniquely by the preceding ones. This 
can be interpreted as an appearance of a singularity where new boundary conditions have 
to be given. Note that this can be neither the initial singularity (which would appear at 
the lowest time values where the initial conditions are fixed) nor a turning point where 
the universe recoUapses (at such a point the coefficients at higher times should be given 
uniquely by the initial value problem in such a way that they are decreasing sufficiently 
fast). It could also mean that at this point our internal time ceases to be a good one; the 
actual interpretation depends on which specific model is considered. 

There are a lot of questions which can be asked about the proposed discrete time 
evolution, most of which can be settled only after studying many models in detail. Maybe 
the first question coming to mind if one is used to the usual evolution in quantum theory 
with a fixed background time is whether the evolution is unitary. But note that this is not 
possible in the usual way because there is an initial and possibly a final time which implies 
that the evolution cannot be described by a unitary operator of the form exp{ih~^tH) with 
a time independent (for a closed system) Hamiltonian H and which would be defined for all 
times — oo < t < oo. First, one would have to adapt this to a discrete time in the form of an 
evolution but this would also imply that each state could be evolved back to negative 
n by applying U"^ arbitrarily often. An initial time together with a unitary evolution 
can appear only in the sense (e.g. in a first order formulation, see below) that there is, 
for any two times rii and n2, a unitary operator £'ni,n2 which describes the evolution from 
time Til to time n2. In order not to be of the form allow an initial time, the 

evolution operators have to be time dependent (although the universe as a whole is a closed 
system), i.e. dependent on rii and n2, not only on the difference. This is indeed the case for 
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our discrete equations because the decomposition of the holonomy multiphcation (e.g. the 
factor in Eq. (^) and the action of the volume operator provide such a dependence 
on n. Note that, whereas in the usual quantum theory time remains forever and does 
not participate in dynamics, in a quantized generally covariant system time is one of the 
dynamical degrees of freedom and can be created as well as annihilated. Therefore, even a 
closed system, which a universe is by definition, must have a time dependent evolution in 
this framework. 

If we have a unitary evolution in the above sense we can define a physical inner product 
in the following way: Let s''^^ and s*^^-* be two solutions of the Wheeler-DeWitt equation. 
In general, i.e. if zero lies in the continuous part of the spectrum of the Wheeler-DeWitt 
operator, their inner product in the kinematical Hilbert space is not defined and both are 
not square summable in the ^^-norm. We can, however, define an inner product on the 
solution space of the constraint by 



,(1) ,(2)^ , ■=Y.Y."^^'^ ■ '^'^ 

=0 L 



for a fixed time ^ f — |- If the evolution is unitary the inner product does not depend 
on the value no and defines the physical Hilbert space Tiphys- Recalling that the evolution 
is generated by a constraint and, therefore, related to a gauge transformation, this inner 
product has the natural interpretation as fixing the gauge by choosing a time no and using 
the kinematical inner product to derive the physical one. In other words, the infinite 
volume of the "gauge" group has been divided out. In this interpretation, we can also 
present the physical inner product in terms of a rigging map: 



.(c'")lc''']:=EE^i?-i,.et,.- 

i=0 L 

Another problem is whether the constraint operator has to be self-adjoint. Accord- 
ing to standard quantization procedures, this should be so because it corresponds to a 
real function on phase space. Nevertheless, in quantum gravity this issue concerning the 



Hamiltonian constraint [|T9| is still debated. We can give here at least a symmetric ordering 



of our Wheeler-DeWitt operators in the form 



2 



where denotes the adjoint in each equal time subspace P„ separately (which is the same 
as * for an operator fixing all these subspaces), and we used the shift operators Sj with 
S* = S-j defined by {Sjc)n '■= Cn+j for j > and {Sjc)n '■= Cn+j if n > — j, {Sjc)n '■= 
if n < j for j < 0. The issue of self-adjoint extensions can, however, be discussed only in 
more explicit realizations. 
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For technical reasons it is interesting to find a first order formulation of the evolution for 
which u = 1. Already for Bianchi models, the order u of the dynamical difference equation 
is rather high (for Bianchi I and IX, e.g., there are maximally two factors of holonomies 
hn in the Euclidean constraint leading to maximally five factors of /i„ in the Lorentzian 
constraint and thus to an order u = 10) which is even increased if we go to isotropic models. 
Note that due to the role of point holonomies in the regularization of the constraint (and 
due to our selection of an internal time) the order of the difference equation increases if 
one introduces some degree of isotropy: now also rotated holonomies can contribute to 
the time spin. Thus, isotropic models have a product of maximally five (flat model) or 
six (closed model) in the Euclidean constraint, and 13 or 15 in the Lorentzian constraint 
leading to an evolution equation of the order 26 or 30. There are two possibihties to reduce 
the order (at least formally): In special cases one could try to take a root, i.e. write the 
equation in the form (A'^ — H'^)c = where 

is the difference operator and H can be interpreted as a Hamiltonian. The second pos- 
sibility, which always works, is to formulate a matrix difference equation for the column 
vectors 

(10) 

\ Cn J 

which have less non-vanishing components for small n, spanning the equal time subspace 
Vn- In both general solution of the first order difference equation leads to evolution 

operators 

Eni,n2'-T^Tn ^ T^n2 (H) 

which describe the evolution from time rii to time n2- 

Assuming such a first order formulation we will now discuss the issue of physical ob- 
servables. 



4 Evolving Observables 

When there is a dynamical formulation of some model system, one is interested in the 
evolution and possible values (in quantum theory spectra and expectation values) of ob- 
servables. Most interesting in the present context is the volume, because we would like to 
know whether the discrete kinematical spectrum is significant also in the physical Hilbert 
space where states solve the Wheeler-DeWitt equation. In the discussion of homogeneous 



quantum geometry [|15| we ignored the Hamiltonian constraint, but now in a quantum grav- 
itational model all observables have to commute with this constraint which is not the case 
for the kinematical volume operator. One can see this by determining the eigenstates of the 
volume operator which, in the dreibein representation, is given by (V"c)n.L = V{n^ L)cn.L 
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(provided the spin network basis by means of which the transition to the dreibein rep- 
resentation has been performed has been chosen as containing only volume eigenstates, 
which is always possible) where the real function V{n, L) is not known explicitly for the 
Bianchi models but depends only on the labels {n,L). We can thus immediately read off 
the eigenstates c^nL^°^ — ^no,nSLo,L to the eigenvalues V{nf), Lq) (possibly degenerate). This 
shows that the (kinematical) volume eigenstates have vanishing coefficients on all but a 
fixed equal time subspace Vn and cannot be solutions of the discrete evolution equation. 

In the framework of classical minisuperspaces (or other generally covariant systems) 
and their standard quantizations there have been proposed several definitions of what has 
to be understood as physical observables subject to a relational evolution with respect to 
another observable |]20|, |T0[ . We use and recall here the method given in Ref . [|T^ because 
it has been applied there to the Bianchi I model such that we will be able to compare this 
construction in a standard minisuperspace quantization with an analogous one in a discrete 
time context. Given a classical function O on the kinematical phase space of the model (so 
that it is not necessarily a physical observable) which does not depend on the canonical 
conjugate of the selected internal time (in an external time formalism such a dependence 
would be impossible, anyway, but here we have to impose this condition because time and 
also its conjugate are chosen from, or functions of, the usual phase space coordinates), one 
first has to quantize it to an operator O acting on wave functions schematically written as 
ilj{a,(f))- As earlier, a denotes the internal time and the remaining degrees of freedom. 
Because O is not assumed to be an observable, Oi/j is in general not a physical state even 
if ip is so. For simplicity, we now make the following assumptions about the evolution: 
a = e* is the (positive) scale factor and the evolution is given in a first order formulation 
by the Schrodinger like equation ih^ = Hi/j with a self-adjoint Hamiltonian H. Note that 
the evolution equation is written in the time parameter t = log a which is unbounded from 
above and below such that there are no obstructions to a unitary evolution with respect 
to this parameter. 

From the operator O one can construct an "evolving observable" 0{t), a family of 
operators depending on the parameter t |]TD[. However, this should not be confused with 
the operator O in the Heisenberg picture (which just would be an equivalent description of 
the same operator); for now both the operator representing an evolving observable and the 
wave function depend on the internal time: the wave function describes the entire history 
of the system which is changed by the action of an operator (a measurement) depending on 
the time when it acts. To construct the evolving observable 0{t) suppose that a solution 
i/j of the Wheeler-DeWitt equation is given. For C'(to) to be a physical observable it has 
to map the physical state "0 to another physical state 0{tQ)tp{a, 0) which will be achieved 
in the following way: first fix the internal time t = to and apply the given operator O to 
the "equal time wave function" ip{a = 6*^,0) (this is possible because O, assumed not to 
depend on the conjugate of time, does not contain derivatives ^) and evolve the resulting 
function to arbitrary times t, i.e. 

(6(to)^)(a = e\ 0) := exp (^-ih-\t - to)H^ d^{a = e*«, 0) . 
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By construction, the right hand side is a physical state and so 0{t) is an observable. 
Note that kinematical expectation values of the operator O are directly related to the 
physical expectation values of 0{t) only if O commutes with the Hamiltonian H where, of 
course, eigenvalues of the evolving observable are functions of the internal time describing 
a relational evolution of some degree of freedom with respect to a. An example for such 
an observable in the above model is the volume which is just given by multiplication with 
e*; as a parameter-dependent multiple of the identity operator it is, however, not quite 
interesting because it describes the "evolution of the volume with respect to the volume" . 
Nevertheless, it shows that for the scale factor all positive real values are allowed, i.e. there 
is no discrete volume spectrum. 

We are now going to present a similar construction for loop quantum cosmology where 
we are, in particular, interested in the implications of discrete time. Our starting point is 
similar to the one above: we have a kinematical operator O and a solution s of the Wheeler- 
DeWitt equation given in the dreibein representation by the coefficients Sn,L where n is the 
discrete time label. Again we assume the evolution to be given in a first order formulation 
by the evolution operators En^^ri2'-T^ni T^n2 of Eq. (|TT|). The same procedure as above 
is to fix the internal time n = no, apply the kinematical operator, and evolve to arbitrary 
time labels: 

{d{no)s)n,L = {Eno,nO s\n=no)L (12) 

where s\n=no is the restriction of s to the subspace Vng. This construction works for any 
kinematical operator which preserves every subspace analogous to the condition that 
O must not depend on the conjugate of time (e.g., in Bianchi models this is the case for 
any operator not containing the holonomy hi which changes the time label n = ji). 

An immediate consequence of the construction together with our discrete time formu- 
lation is that only discrete values for the parameter n in evolving observables 0{n) are 
allowed, so evolution is always with respect to a discrete time. Most interesting in cosmo- 
logical models is the volume which is kinematically quantized by the operator V and has 



a discrete spectrum [|T5|. We can apply the construction of an evolving volume operator 
V{n) which describes the evolution of the volume with respect to discrete time n (because 
now time is not chosen to be the scale factor, the relational evolution V{n) is non-trivial). 
Most important, also for the full theory, is the question of whether the spectrum of this 
evolving volume is again discrete and somehow related to the kinematical spectrum. 

In order to make the relation between kinematical and physical spectra more precise, 
we now introduce some formalism. We write the evolution as an operator En., which, acting 
on a state d in the "initial" subspace P„, yields a complete history given by a physical 
state: 

Eno,- - "^no ^ "^phys ) (-^no,-^)"i^ (-^no,n'^)L • 

Because {Eno,.d)\n=no = d in general and £'„y_.(s|„=„o) = s for a physical state s, En^,- is 
the inverse of the projection operator 

TTriQ • "^phys T-^uq ; S I > S | n=no 
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on physical states in TYphys- On the kinematical Hilbert space T^km? is not 
invertible but we can easily calculate its adjoint 

which is just the inclusion map of as a subspace of Tikin- From the equation 

n,L L 

we directly infer Tr*^^ = (both operators are bounded). 

We can now write an evolving observable constructed from a kinematical operator O 

as 

d{n) = En, o 6 o 7r„ . (13) 

Of physical interest are the expectation values (s, C'(?T,)s)phys which, provided s is a physical 
state, describe the relational evolution of the observable O with respect to discrete time n 
in the given history. These expectation values are directly related to the physical spectrum 
of O. On the contrary, the kinematical spectrum of O is directly related to the expectation 
values {d, Od) of the kinematical operator O in an equal time subspace P„. Choosing a 
physical state s with 7r„s = d and using the above definitions and adjointness relations, we 
can write this as 

{d,dnnS) = {d,'KnEn,d'KnS) = {L„d,d{n)s) 

and see that the kinematical expectation values are given by matrix elements of 0{n) with 
respect to an unphysical state Lnd and a physical state s. This immediately shows that the 
kinematical expectation values do not have any physical meaning. 

Although individual expectation values in a given physical state are unrelated to kine- 
matical ones and a computation of an evolving expectation value (s, 0(n)s)phys in a given 
history s can be very complicated, it is simpler to study the set of possible outcomes of 
expectation values, i.e. the spectrum, of an observable 0(no) at a fixed parameter uq. To 
that end, we proceed in the following way, using the example of the volume operator: We 
fix an no and diagonalize the kinematical volume operator V on the equal time subspace 
VnQ. Because we need a first order formulation in order to define V^uq), we assume 
to be the space of w-column vectors of Eq. (|l^). Eigenstates of V in are denoted 
as < i < uj — 1, and fulfill the eigenvalue equation Vv^'^'^^ = V{nQ — |, L)v^^'^'^ with 

the kinematical eigenvalues V{n,L). 

We can then evolve each of these eigenvectors to complete histories s('^o,2,i) -which form 
a complete set in the physical Hilbert space. By definition, we have 

y(no)s("«'^'^) = En,,Vv^''''^ = En,,V{n, - f , L)t;(*'^) = V{n, - f , Ds^'^"'*'^) (14) 

demonstrating that the kinematical eigenvalues are also the physical ones. Of course, this 
does not tell anything about evolution because this is coded in the dependence on n of 
expectation values of V{n) in a fixed physical state (the set of eigenstates used above is 
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n-dependent). But it shows that kinematical spectra are relevant for physical operators. 
We have here the first models where a discrete spectrum of a metrical operator in the 
physical Hilbert space has been derived. Note, however, that our conclusion depends on 
the assumption that the spin n = ji, or another combination of the spin labels, defines a 
sensible time variable. 



5 Discussion 

Although we did not discuss any particular model but concentrated on conceptual issues, 
it has by now become clear that our proposal of loop quantum cosmology is very different 
from the standard minisuperspace quantum cosmology. In view of the discrete structure of 
space revealed in loop quantum gravity this is what one would intuitively expect because 
then also time should be discrete implying departures from the conventional continuous 
time of quantum mechanics. We emphasize here that, whereas the standard quantiza- 
tions are equivalent to the treatment of (relativistic) quantum mechanical systems and 
heavily rely on methods which are assumed to be inapplicable in a full quantization of 
General Relativity, our attitude was always to be as close to the general framework of 
the loop quantization of General Relativity as possible. In this sense, we regard results 
of loop quantum cosmology as being more trustworthy compared to an extrapolation of 
minisuperspace results to the full theory. Manifestations of the close relationship to full 
loop quantum gravity are the discrete geometric spectra and the very similar Hamiltonian 
constraint operators. This lead us to the derivation of a discrete time evolution and of dis- 
crete physical spectra of geometric operators as new contributions to the loop quantization 
program. 

In contrast to standard minisuperspace quantizations we did not perform the symmetry 
reduction directly at the classical level, but were able to interpret the states as symmetric 
states in the kinematical sector of a quantization of the full theory leading to a preservation 
of discrete geometric spectra. But still, the symmetry reduction is very restrictive which 
can be observed in the phenomenon of level splitting: whereas the full volume spectrum is 
almost continuous for large eigenvalues (they are given by arbitrary sums of irrational num- 
bers, the distance between subsequent ones can be made arbitrarily small at sufficiently 
high values), the reduction to homogeneity leads to a degeneration probably undoing the al- 
most continuity (the spectrum is not known explicitly, however); for the highest symmetry, 
isotropy, the distance between subsequent levels even increases with increasing eigenvalues 
p!5| . Nevertheless, the relative difference decreases as + 1)2 —j^) ~ which still 



can be sufficient to recover the usual classical continuous space for large volumes (com- 
pare with energy levels of the harmonic oscillator). Anyway, quantum cosmological models 
are most interesting for small volumes, i.e. close to the classical singularity, where the 
discreteness of both the full and the symmetric spectra can be expected to be relevant. 
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5.1 Comparison with the Full Theory 

Of course, as compared to the full theory we are in a very special situation if we are 
studying the dynamics of homogeneous models because, first and foremost, we have a 
natural candidate for an explicit internal time. This fact enabled our discussion of the 
dynamics of loop quantum cosmology. Another simplification comes from the reduction of 
degrees of freedom leading to models which are defined on a fixed simple graph. Future 
analytical and numerical computations will benefit from this reduction. 

On the other hand, Wheeler-DeWitt operators of homogeneous models P] are very 
similar to a single vertex contribution of the operator in the full theory 0. Therefore, 
their analysis will not be much simpler. However, even without having explicit solutions 
of the Wheeler-DeWitt equation we were able to derive consequences for the evolution, 
which hopefully will help to shed light on the dynamics of the full theory. 

5.2 Difference vs. Differential Evolution Equations 

The main discrepancy between standard minisuperspace quantizations and the one pre- 
sented here is that time evolution is described by a differential equation for the former and 
by a difference equation for the latter ones. For a numerical analysis this makes no dif- 
ference because differential equations, if they cannot be solved analytically, are discretized 
anyway. 

Conceptually, the following question may be more interesting, namely whether our 
discrete time evolution equations can be interpreted as a discretization of a continuous 
(semi classical) time evolution. If this is possible, the continuous time formulation would 
be of relevance for a discussion of the classical limit. However, this problem is highly non- 
unique which implies that there may be different continuum pictures described by one and 
the same discrete evolution equation. 

Such a behavior is also suggested by the problem of consistency which often has to be 
faced in numerics (see e.g. Ref. 0]). If a differential equation is discretized for a numerical 
analysis, it is possible that a solution converges to a solution of another differential equation 
if the discretization is refined. In particular for higher order approximation schemes, there 
are additional and often very weird solutions which have to be under control in a sensible 
code. 

Recall that our discrete equations (H) are usually of a very high order as compared to 
the order two of the differential Wheeler-DeWitt equations. Although by itself this is no 
obstruction to recover a standard low order differential equation, it implies that we have to 
expect in general very many solutions which can have no classical counterparts. Note that, 
contrary to the numerical analysis of a differential equation where additional solutions of 
its discretization have to be suppressed, in our context the discrete description is regarded 
as being more fundamental. Therefore, additional solutions have to be suppressed only 
in semiclassical regimes; they are genuinely quantum ones and may be of relevance for a 
discussion of quantum modifications of the classical singularity. 
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5.3 Avoiding the Classical Singularity 



The high order of our discrete equations has a consequence for the question of whether 
there is always a positive probabihty for a physical history to have a singularity. Under 
"singularity" we will in this context understand a geometry with vanishing volume, i.e. 
we ask whether there is for each physical state s an uq such that V{nQ)s = 0. Because 
the evolution is linear, this is the case if the kinematical volume operator V applied to an 
equal time slice of the physical state s is zero: VTTnoS = with 7r„gS ^ 0. 

On the classical side, there are the singularity theorems |^ which state that any 
classical space-time has a singularity under some conditions on the matter content (and 
the cosmological constant). In Bianchi models a singularity is generic, whereas isotropic 
models have a singularity only if matter is coupled which fulfills a certain positive energy 
condition (which may also have to compensate a positive cosmological constant). Without 
matter, spatially isotropic space-times are maximally symmetric (Minkowski, DeSitter or 
Anti-DeSitter space) and thus cannot have a singularity (except for coordinate singularities 
which are introduced by an inappropriate slicing or time coordinate). The question is now 
whether there is a quantum analog of these theorems or whether a singularity can be 
avoided there in the sense that the probability for a universe to be in a singular geometry 
vanishes. 

In Bianchi models, the kernel of the volume operator is very large, even for positive 
times n > such that a discussion here is impossible without explicit solutions of the 
Wheeler-DeWitt equation in a specific model. 

In isotropic models, however, there are only three states with zero volume ||TB[. The 



order of the discrete evolution equation, on the other hand, is much higher so that we can 
easily demand that the amplitude Sn,L, where L denotes the insertion and matter labels, of 
a physical state vanishes in these degenerate states thus avoiding the singularity. Note that 
this is independent of the matter coupled to the isotropic model as long as time is made 
only of gravitational degrees of freedom as we did above. Thus we see, that also for matter 
which classically inevitably leads to a singularity a degenerate geometry can always be 
avoided in quantum solutions of isotropic models. However, as already discussed, the high 
order discrete equations governing the quantum time evolution have a wealth of solutions 
most of which will not correspond to classical ones. Therefore, also in isotropic models 
an avoidance of the singularity can be concluded only if one finds a solution which has 
vanishing amplitudes on the degenerate geometries and corresponds to a classical solution 
for large volumes. 
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